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ABSTRACT 

Soft X-ray absorption in excess of Galactic is observed in the afterglows of most gamma-ray bursts (GRBs), 
but the correct solution to its origin has not been arrived at after more than a decade of work, preventing its 
use as a powerful diagnostic tool. We resolve this long-standing problem and find that absorption by He in the 
GRB's host H II region is responsible for most of the absorption. We show that the X-ray absorbing column 
density (Nh x ) is correlated with both the neutral gas column density and with the optical afterglow's dust 
extinction (Ay), This correlation explains the connection between dark bursts and bursts with high Nh x values. 
From these correlations we exclude an origin of the X-ray absorption which is not related to the host galaxy, 
i.e. the intergalactic medium or intervening absorbers are not responsible. We find that the correlation with the 
dust column has a strong redshift evolution, whereas the correlation with the neutral gas does not. From this 
we conclude that the column density of the X-ray absorption is correlated with the total gas column density 
in the host galaxy rather than the metal column density, in spite of the fact that X-ray absorption is typically 
dominated by metals. The strong redshift evolution of Njj x /Ay is thus a reflection of the cosmic metallicity 
evolution of star-forming galaxies and we find it to be consistent with measurements of the redshift evolution of 
metallicities for GRB host galaxies. We conclude that the absorption of X-rays in GRB afterglows is caused by 
He in the H II region hosting the GRB. While dust is destroyed and metals are stripped of all of their electrons 
by the GRB to great distances, the abundance of He saturates the He-ionising UV continuum much closer to 
the GRB, allowing it to remain in the neutral or singly-ionised state. Helium X-ray absorption explains the 
correlation with total gas, the lack of strong evolution with redshift as well as the absence of dust, metal or 
hydrogen absorption features in the optical-UV spectra. 

Subject headings: gamma-ray burst: general — early universe — dark ages, reionization, first stars — galaxies: 
ISM 



1. INTRODUCTION 

Long-duration gamma-ray bursts (GRBs) accompany the 
death s of some stripped-e nvelope massive stars (Hjor th et alJ 
12003b IStanek et alll2003l) . And while they have been ac- 
knowledged as excellent probes of both the typical interstel- 
lar medium of their host g alaxies and foreground absorbers 
(e.g. Jakobsson et al. 2004), their use as probes of the progen- 
itor star's circumstellar material has met with limited success. 
The main difficulty in using GRBs to examine their immediate 
environs is that they are so luminous across all wavelengths 
that they can be expected to ionis e gas and destroy dust to 
very large distances from the bur st (Waxman & Drainel l2000[ 
iFruchter et "ail 120011: iPerna et all 120031) . potentially eradicat- 
ing all traces of the circumstellar environment from their spec- 
tra, except potentially in hot gas. The best hope for observing 
the environment near the burst could therefore be at X-ray 
wavelengths. 

From the early soft X-ray observations of GRBs, large 
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absorptions were observed, initially believed to be con- 
sistent with the light from the afterglow making its way 
out of a molecular cloud and su ffering absorption from 
the metals in the molecular cloud (IGalama & Wiiersl 120011: 
iReichart & Pried 12002b ICampana et al.ll2006al) . However, it 
has become clear since then that the picture is not that sim- 
ple. Indeed, the distribution of X-ray absorbing column den- 
sities is substantially larger than one would expect in such 
a scenario. And it has been shown that the X-ray absorp- 
tion is not readily connected to either the neutral hydrogen 
column observed in the optical whether corrected for metal- 
licity or not, to low-ionisation gas, or to the column of dust 
( Wats on et al.l2007t ICampana etaLll2010b iSchadv et a0201 lh 
Zafa r et al.ll201llh 

As we progress in our understanding, the nature of the X- 
ray downturn at soft X-rays seems to become less and less 
clear. We still do not know how it arises, and indeed we are 
beginning t o question whether it is even due to photoelectric 
absorption (Butler & Kocevski 200 7]) . This question has been 
a deca de in the making (IGalama & Wiiers 2001 ; Watso net al.1 
2002) and is one of the outstanding observational issues re- 
lated to GRBs. Many suggestions have been made as to its 
origin, among them highly-ionised gas in the intergalactic 
medium (IGM), the host galaxy molecular cloud, or circum- 
stellar material related to the GRB progenitor. 

If we can resolve the origin of the X-ray absorption, we may 
find out about the immediate environment of the GRB, and 
hence know about its progenitor and the matrix in which it 
was c reated; or, if recent suggestions are correct (Beh ar et alJ 
1201 lh . resolve the missing baryons problem by detecting the 
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warm-hot intergalactic medium. 

In § |2] we present the problematic interpretation of the ex- 
isting data and discuss how the X-ray column density is de- 
termined and presented. In § [3] we present our results on the 
correspondence between the X-ray absorbing column density 
and the gas and dust columns. The origin of the X-ray ab- 
sorption is explored in § |4] while its specific properties are 
discussed in § [5] The implications of our findings are exam- 
ined in §|6]and our conclusions presented in §|7] 

2. OBSERVATIONAL FACTS ABOUT X-RAY ABSORBERS IN GRBS 

We now array the facts we know about the X-ray absorption 
in GRB afterglows. 

In the Galaxy, there are strong correlations between 
the neutral hydrogen column density, Nn, the dust 
column density, A y, and the X-r ay abs orbing col- 
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iRachford et all 120091: iGiiver & Ozeil 120091 iWatsonl 
201 ll) . The proportionality constant is different in the Mag- 
ellanic Clouds for Nn/Nu x , but corrected for the lower 
metallicities of these galaxies, the ratio is quite simil ar 
dKoornneefj|1982t iBouchet et alj[l985l i Fitzpatricld 1 1 985allbL 
1986b [Martin et al lll989t iGordon etaLll2003t iDobashi et al] 
20081 IWeltv et all 120 121) . The metals-to-dust ratio, N Hx /A v , 
is relatively const ant not only in local group galaxies 
(IWeltv et alJ 120121). but also in galaxies at cosmologi - 
cal redshifts (e.g. iDai et all l2006t iDai & Kochanekl I20091) . 
However in GRB afterglows, the comparison with H I 
and Ay shows that A^h x is significantly larger than ex- 
pected, whether metallicity corrections are included or not 
(IZafar et alJ 12011b ICampana etaD 120121) . N Hx /A v is typi- 
cally an order of magnitude larger in GRBs than observed 
in other galaxies with a wide variation (IWatson et al.l 120061 
12007b iSchadv etaLlHoTob IZafar et all 120 111) . Finallv. it was 
recently shown that the Nu^/Ay ra tio evolves strongly with 
redshift (IWatson & Ja kobsson 120121) . This is strange and so 
far remains unexplained. We explore this phenomenon below. 

The consensus view has so far been that the bulk of X-ray 
absorption is not due to what is causing the dust extinction 
and gas absorption observed in the opti cal and is probably 
therefore due to very h i ghly ionised gas (IWatson et al.H2007b 
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2.1. Practical considerations 

The dow nturn observed in the soft X-ray spectra of GRB 
afterglows dGalama & Wiiersll2001b IWatson et all 120021) was 
initially assumed to be the same as that observed in the Galaxy 
and in most extragalactic sources: photoelectric absorption 
due to metals along the line of sight. This absorption is pri- 
marily due to the inner shell electrons of the most abundant 
metals. In the energy range observed by most modern X-ray 
detectors (approximately 0.2-10keV) these metals are in par- 
ticular C, N, O, Ne, an d S with cont ributions from the L-shell 
of Fe (see for example [Wilms et"al]l2000l) and a contribution 
from He. Oxygen provides approximately 40% of the opacity 
in the Swift-XRT's, passband at z = 0. Because the absorption 
is caused by the inner shell electrons, mildly ionised gas still 
absorbs soft X-rays, as do atoms in dust or molecules, with 
relatively little alteration. Thus, soft X-ray absorption is a 
means to determine the total column density in metals in front 
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Figure 1. Soft X-ray spectrum of the afterglow of GRB 121027A, convolved 
with the XMM-Newton instrument response. The best-fit models for a power- 
law with a neutral, solar metallicity absorber (black, solid line), and a pure 
He absorber (red, dashed line) both at z = 1.773 are shown. The models are 
practically indistinguishable. 



of an object. Unfortunately, the low effective spectral resolu- 
tion to absorption features of all of our soft X-ray spectra is 
such that individual features cannot be distinguished (includ- 
ing the spectra obtained with XMM-Newton and Chandra's 
gratings, since the total numbe rs of counts obtained with these 
instruments is so low, e.g. jMirabal et aD 12003b [Butler et al. 
1200112001: ICampana et al.ll201 lh 

It's worth noting here that at low metallicities, the absorp- 
tion is increasingly dominated by He (and to a lesser extent, 
H, if it exists in the neutral phase), and there is no clear way 
to distinguish absorption by different elements. As an illustra- 
tion, in Fig.Q] we show the s pectrum of the X-ray afterglow of 
GRB 1 21027 A at z= 1.773 ilTanvir et al.ll20H iKruhlereFall 
l2012al) . fitted with a traditional solar abundance absorber, and 
with a pure He absorber. The difference between the fits for 
the models is negligible. However, the total column density 
in He is nearly a factor of seven larger for the pure He ab- 
sorber. For singly-ionised He, this is another factor of ~ 25% 
larger again, since the cross-section for neutral He is some- 
what larger in soft X-rays than for He + . This means that for a 
H II region, with a considerable fraction of He II but no met- 
als, the gas column density would be about 8 times the column 
required for a neutral, solar metallicity gas at this redshift. 

The lack of spectral resolution good enough to discrimi- 
nate individual features gives rise to complications and con- 
fusion in the presentation of results. The first is what column 
density to report. Since we typically not only do not know 
the state of the matter or indeed the precise elements we are 
observing, a column density of any given species cannot be 
reported. For this reason an equivalent column density in hy- 
drogen is usually reported assuming, typically, a 'solar' abun- 
dance of the elements — usually the default abundan ce set in 
the software Xspec, from An ders & Grevesseldl989l) . In spite 
of newer estimates of the solar abundances that are about 40% 
lower in metals (e.g. lAsplund et al.l 120091) . it seems that the 
lAnders & Grevesseld 19891) met allicity is a better approxima- 
tion than Asplund et al. (2009) for a typical Galactic sight 



line (IWatsonl 20111) . For this reason and for ease of compar- 
iso n with previous estimates, it makes sense to continue to 
use lAnders & Grevessd (Tl989) metallicities unless we know 
what the absorber is. The second complication is that the ab- 



3 



sorbing column density is occasionally reported as if the ab- 
sorber was at z = 0, typically where the redshift is unknown. 
The absence of spectral features means that the redshift must 
be determined from observations at another wavelength, typi- 
cally optical or UV. Since the absorption is shifted out of the 
bandpass with redshift, the observed column density drops 
with redshift approximately as (1 +z)~ 26 . Naturally the col- 
umn density reported at z = is thus substantially smaller than 
the actual value determined at the correct redshift. Third, the 
value found for the column density is almost always deter- 
mined assuming a neutral column density. Ionised gas has a 
smaller absorbing cross-section, and indeed, a different ab- 
sorbing pattern. 

3. HIGH X-RAY ABSORPTION AT HIGH REDSHIFT: WHY DOES 
N Hx /A v EVOLVE WITH REDSHIFT? 

It was noted by ICampana et al.l (120101) that there was a 
correlation between the X-ray absorbing column densities 
in GRB afterglows and their redshift, i .e. the highest red- 
shift bursts were the most absorbed. iBehar et ail (1201 lb 
demonstrated convincingly that this was not an artifact of 
the fitting and assumptions, and suggested that this relation 
was due to increasing foreground absorption by the highly- 
ionise d metals component of the i ntergalactic medium. How- 
ever, IWatson & Jakobssonl (120 1 2|) found the ratio of the X- 
ray absorption to dust extinction Nn x /Ay in GRB afterglows 
evolved with redshift, and that it was this evolution that an- 
swered the mystery of the dearth of high X-ray absorbed 
GRBs at low redshift. This answer was in many ways 
more mysterious than the question. Why would the X-ray- 
determined metal column density increase relative to the dust 
column density with redshift? Here we examine the relation- 
ship between Nh x and Ay as a function of redshift to under- 
stand the origin of this peculi ar evolution using (a slightly 
enlarged version of) the data in lWatson & Jakobssonl (120 121) . 
These data are presented in detail in Wa tson et al. (in prep.). 

A c lose examination of Fig. 3 in Wa tson & Jako bsson 
(120121) hints at a possible correlation between the metal col- 
umn densities and the extinction that evolved with redshift. 
While there is no immediate correlation apparent in the full 
dataset, splitting the data by redshift clarifies the situation. 
In Fig. |2j the distribution of X-ray absorbing column density 
with extinction is plotted for GRBs in four redshift intervals. 
The correlation then becomes more apparent, with correla- 
tion probabilities of greater than 99.9%, 97.5%, and 96% for 
GRBs in the intervals z < 1, 1 < z < 2, and 2 < z < 4 re- 
spectively based on Kendall's r. The combined probability 
that such correlations occur randomly is less than 1 x 10~ 5 . 
Similar probabilities are obtained using Spearman's p. While 
the correlation is highly significant, there is clearly significant 
scatter in this correlation. In sum: the X-ray absorption is cor- 
related with the dust extinction and the relation evolves with 
redshift. 

3.1. Origin of the correlation between Nn x and Ay 

At first glance, the fact that dark bursts are preferen- 
tially associate d with high jV H y G RBs dKriihleret alj|201 lb 
ICampana etaDl2012b iFynbo et al.1 120091) implies an N Hx -A v 
correlation. The evolution with redshift and the sc atter how- 
ever, was concealing it ("e.g. ICampana et a"fll2010t) . The first 
interpretation might be that this is simply host ISM in both 
dust and metals, resulting in high Nn x an d Ay values. How- 
ever, we know that the absorption in the optical-UV is dom- 
inated by low-ionisation gas and dust at distances of hun- 
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Schadv et all 1201 It IVreeswiik et al.l 
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that these low-ionisation metal column densities correlate 
well with the observed Ay (Zafar et al. in prep., see also 
Ide Ugarte Postigo"et al. 2012). For the observed Nn x -Ay cor- 
relation to hold, the cool, ISM would have to dominate the 
colum n density. And indeed, it has been argued ( Kriihl eret al.1 
|20TTI) that even with both components entirely uncorrected, 
an approximate correlation between H I and Nu x would be 
found at high column densities if the H I column is sufficiently 
dominant. 

However this cannot be the origin of the correlations we 
observe with the X-ray column. First, the H I column den- 
sity, corrected for metallicity, is almost never as large as the 
observed X-ray column d ensity (see Fig. [3] also Zafar et al. 
1201 ltlKrtihler et aT1l201 ll) . and the correlation extends to low 
column densities. We show in Fig.|2]the X-ray abs orption val - 
ues after subtraction of 2.0 x 10 21 cm" 2 x Ay (IWatsonll20lTb . 
which is the absorption expected for the metals typically as- 
sociated with such a column of dust (Zafar et al. in prep.). 
The remaining X-ray absorption is still correlated with Ay, 
with probabilities of 96%, 94% and 92% for the three lower 
redshift sets. The combined probability of these datasets all 
having such strong correlations by random chance is only 
- 3 x 10~ 4 . 

Second, the correlation between the hydrogen column den- 
sity and the X-ray absorbing column density i s tighter before 
accounting for metallicity (see Figs. 3 and 4 in Campa na et al.l 
2010). This shows that the correlation of the X-ray absorbing 
column density is primarily with the full gas column density 
and not particularly with the neutral metals. Third, there is 
no obvious reason there should be an evolution in the Nn x /Ay 
ratio with redshift in this scenario. We conclude, therefore 
that the typical host galaxy ISM is not responsible for the cor- 
relation and that there is a correlation between the dominant 
component of the X-ray absorbing column density and the to- 
tal gas column density. 

3.2. High X-ray column density, high gas column density 

Such a correlation is odd at first, since the X-ray absorption 
is usually principally due to metals. However, what it is telling 
us, is not that the gas column is responsible for absorbing the 
X-rays, but that where the X-ray column density is high, so 
too, often, is the gas column density. This correlation with the 
gas column density then explains the strong redshift evolution 
of the Nn x -Ay relation, but the apparent lack of redshift evolu- 
tion of the Nh x -Nhi relation (Pig. [3j. It is simply the average 
metallicity evolving with redshift. At high redshifts we have 
high X-ray column density objects with high gas column den- 
sities. These objects are significantly metal poorer than at low 
redshift, and hence (via the virtually constant dust-to-metals 
ratio demonstrated in Zafar et al. in prep.), have a lower Ay. 
In other words, the evolving Nn x -Ay relation is a reflection of 
the cosmic metallicity evolution. 

We can now draw a few important conclusions from this. 
First, the X-ray absorption component must be associated 
with the host galaxy since it correlates with the gas column 
density in the host. Thus it cannot be due to the intergalactic 
medium, for example. Second, since the maximum values of 
7Yh x are very similar at high and low redshift, and the mean 
distribution certainly does not decrease to high redshift, as 
would be expected if we were observing an X-ray absorbing 
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Figure 2. X-ray absorption versus extinction in GRB afterglows. The four panels present GRBs in different redshift intervals: top left z > 4, top right 2 < z < 4, 
bottom left 1 < z < 2, bottom right z < 1 . The solid line is the best-fit, fixed-slope line to the data. The dashed lines mark the approximate limits of the 
metals-to-dust ratios reported for the local group of galaxies. The metal and dust columns are clearly correlated in each redshift interval (P > 95% in the three 
z < 4 datasets). The best-fit metals-to-dust ratio rises at high redshift: from ~ 1 X 10 22 cm -2 mag -1 Ay at the lowest redshift to ~ 6 X 10 22 cm -2 mag~' Ay at the 
highest. The datapoints in faded colours are Nn x —(2.0 X 10 2 ' cm -2 X Ay) versus Ay. The remaining X-ray absorption when the metal absorption corresponding 
to the near-universal metals-to-dust ratio (Zafar et al„ in prep.) is subtracted (2.0 X 10 2 ' cirT 2 X Ay) is still clearly correlated with Ay, showing that the Nn x -Ay 
correlation we observe is not due to the neutral medium metals associated with the dust forcing a correlation at high Ay values. 



column dominated by metals in the ISM gas, it suggests that 
we are either observing metals ejected by the progenitor star 
itself or that a component of the gas largely unaffected by cos- 
mic metallicity evolution is the absorber. 

3.3. The explanation of the ' dark burst' -N^ connection 

The overall correlation between the X-ray column density 
and the dust extinction can be explained by a model where 
GRB progenitors residing in the hearts of galaxies are sur- 
rounded by higher density ISM. Both H II regions and WR- 
nebulae are known to expand to smaller sizes in the hearts 
of galaxies than those on the outskirts o r outside the plane 
dHunt & H irashita 2009; Stock & Barlow 2010). These more 
compact sizes result in higher observed column densities 
through the nebula or H II region. GRBs occurring in high 
ISM density regions will therefore have higher X-ray column 
densities. Naturally, the afterglows of such bursts are also 
much more likely to encounter high density ISM sightlines on 
their way out of their hosts. The approximate correlation ob- 
served between the X-ray column density and the dust extinc- 



tion Fig.|2]is then due to this effect. The fact tha t 'dark bursts' 
have on average highe r X-ra y column den sities (iKriihler et al. | 
201 U [ Camnana et al.ll20Tl iFynbo et"aT]l2009tlKruhler et al.1 



2012b) is clearly an observational corollary. 



We can confirm this explanation if we observe an overall 
and rough correlation between the X-ray absorbing column 
density and the neutral hydrogen column density. This is the 
case. In Fig. [3] we plot the X-ray absorption against the neu- 
tral hydrogen column density determined from the Lya line. 
There is an overall correlation between the two (96% confi- 
dence including all the data), with a few strongly discrepant 
outlying objects with very low H I column densities — which 
in this scenario are objects that happen to find a clean line of 
sight out of the host. In an earlier, smaller sample, we found 
no evidence for a correlation between the X-ray and H Lya 
column densities, pa rtly due to these few very low H Lya col- 
umn density objects (IWatson et al.ll2007h . 

The mysterious dr op in the observed dust col umn for higher 
redshift sources (see Wat son & Jakobssonl2012l) is readily ex- 
plained in this context. If we assume that the average metallic- 
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Figure 3. Gas (H I) column density as a function of X-ray absorbing col- 
umn density. Objects at high (z > 3) and low (z < 3) redshifts are marked 
with squares (blue) and circles (black) respectively. At the higher values of 
A^hi, i.e. for most objects, the gas and X-ray column densities are correlated, 
though with a considerable intrinsic scatter. The correlation is similar to that 
observed between the dust extinction and the X-ray absorption (Fig. [2). How- 
ever, unlike the dust correlation, there is no obvious evolution with redshift. 
Note that due to the low energy cut-off, the detectability threshold for the 
X-ray absorption is strongly redshift dependent. This detectability effect is 
clear with the high redshift objects (black) having higher X-ray absorption 
and upper limits. 

ity of galaxies is low at high redshifts and the gas columns are 
approximately similar for similar X-ray absorptions at all red- 
shifts, and the dust-to-metals ratio is roughly constant, then 
the apparent increase in the Nn x /Av is simply a result of the 
decreasing metallicities at these high redshifts. If this expla- 
nation is correct, we should observe an evolution of increasing 
Nh x /Av with redshift, but no such evolution in the Nm/Nn x 
ratio. Again, this is what we observe, with no apparent evolu- 
tion of the N m /N Hx ratio (Fig.E). 

A potentially fascinating tool resulting from this discov- 
ery is that the observed evolution of the mean Av/Nh x for 
GRBs can therefore yield an approximate indication of metal- 
enrichment of star-forming galaxies with time. At first glance, 
for example, we can argue that the metallicity in star-forming 
galaxies decreases by approximately a factor of 3 between z < 
1 and 2 < z < 4 and by a factor of 5 between z < 1 and z > 4. 
This is a flatter sl ope to the metallicity evolution than found 
with QSO-DLAs dProchaska et al.1 120031) . GRB host metal- 
licities, on the other hand, where they have been measured, 
also show an evolution with redshift, but the evolution is sig- 
nificantly slower than for QSO-DLAs and is quite consistent 
with what we observe here (iFynbo et al.l2006tlLeyesaue et al.1 
I2010tlde Ugarte Postigo et al.ll2012tlThone et al.ll2012l) . This 
is reassuring for the interpretation presented above. 

4. WHAT IS THE X-RAY ABSORBER? 



Given the correlation with gas in the host galaxy, it is clear 
that the X-ray absorber must have some connection to the host 
ISM, but not the typical free, low-ionisation ISM observed far 
from the burst. We are thus left with two possibilities: the 
immediate cloud in which the burst progenitor was born, or 
circumstellar material. 

4.1. The host molecular cloud 

It has been proposed that the X-ray absorption is 
related to gas and dust in the molecular cloud in 
which the GRB progenitor is formed, but that the GRB 
ionises the gas and destroys the du s t to potentially large 
distances dWaxman & Draind 120001: iFruchteret all 120011 
iGalama & Wiiersl I200U iPerna et all 12003b . leaving the X- 
ray absorption by highly-ionised metals as the only trace. 
However two pieces of evidence argue against this possibil- 
ity. First, the low-ionisation optical absorber is observed to 
lie at very large distance s in almost every case , hardly ever 
at dist an ces of a few pc (IVreeswiik et a l. 2007; D'Eli a etldl 
l2009allR l20Tlt iLedoux et al.ll2009i: iFox et alj|2008l). Molec- 
ular hydrogen is very rarely detected (Ledoux et al.l 120091 
iProchaskaetaLl 120091 IWhalen et al.1 12008|), thoughthere is 
a strong dust-extinction bias against such detections. And 
there is no one-to-one correlation between the X-ray absorp- 
tion a nd the neutral hydrogen column density (IWatson et al.1 
120071) . This means that the hydrogen related to the X-ray ab- 
sorber must be dissociated and ionised. However there are 
not enough UV continuum photons in most GRBs to dis- 
sociate and then ionise the hydrogen out to distances much 
greater than ~ 1 pc from the burst for typical inferred col - 
umn densities (Wax man & Draing|2000t IWatson et al.l 2007). 
Second, there is little evidence for the moderately ionised 
species one might expect to exist, lying in the zone between 
the highly-ionis ed X-ray absorber and the almost neutral op- 
tical absorber (ISchady et al.ll201 ll) . Ultimately, the GRB is 
not UV-luminous enough to ionise the hydrogen to super-pc 
distances, while the molecular clouds in which massive stars 
form should often be larger than this size, and we do not ob- 
serve any of the signatures of such molecular clouds, of their 
ionisation by the GRB, or of the remnant neutral/molecular 
material we would expect if th e GRB were 'burning' its way 
out of its molecular cloud (see Wa xman & D raine 200(|). We 
therefore conclude that host molecular clouds cannot be the 
origin of the X-ray absorption in the general case. 

4.2. The host H II region 

We must be careful to distinguish between a molecular 
cloud and what once was a molecular cloud, which, after only 
a few million years will be substantially ionised by the mas- 
sive stars formed in it, transforming it into a H II region, likely 
before the burst occurs. We are then discussing H II regions 
rather than molecular clouds, and the arguments related to 
ionisation of hydrogen by the GRB do not apply. The to- 
tal metal column densities are in the right range to explain 
the observed X-ray absorption, in H II regions with sizes up 
to several tens of pc in size (Fig. 4. iHunt & H irashita 2009). 
However, the principle objection to this scenario is the obser- 
vation that the A^h x Mv ra ti° evolves with redshift, which can- 
not readily be explained in this scenario unless the absorption 
is not due to the metals in the gas. 

The absence of an observed decrease in the X-ray column 
density to high redshift is a strong indication that the X-ray 
absorber is not related to the metals in the general ISM of 
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the host galaxy. If it was, the overall evolution of the cosmic 
metallicity would imply a decrease of two orders of magni- 
tude in the absorption to z = 8. If the X-ray absorption was 
influenced by the mean ISM metallicity, we would anticipate 
a strong decrease in the mean X-ray absorption. If anything, 
the mean X-ray absorption appears to increase w ith redshift 
dCampana et all 120121: IWatson & Jakobssonl 120121) and this 
cannot be a selection bias, since the total number of objects 
at high redshift would become unfeasibly large if there was 
a missing population of low-absorption systems at high red- 
shift. In addition, there should be no strong bias against dis- 
covering (not measuring) low-absorption systems at high red- 
shift. We cannot see how an absorber dominated by the mean 
metallicity of the host galaxy can reproduce the observed high 
absorption systems at high redshift, and therefore are com- 
pelled to exclude metals in the host H II region as the ori- 
gin of the X-ray absorption. However, it is worth considering 
whether the He in the host H II region could provide enough 
X-ray absorption. In that case we would need to demonstrate 
that H II regions have the correct distribution of gas and dust 
column densities, and radii, and we would have to explain the 
lack of observations of moderately-ionised species. 

4.3. Ejecta from the progenitor 

It is natural enou gh, considering the association of GRBs 
with type Ic SNe dGalama et al] [19981: IStanek et all 12001 
Hiorth et al. 2003), to suggest that GRBs have heavily mass- 
losing stars (specifically, Wolf-Rayet (WR) star s) as progen- 
itors (e.g. lWooslevl[l993t iChevalier &"Dlll999l) . Recent ob- 
servations of WR stars in the Galaxy and Magellanic clouds 
show that the WR nebula, filled with the highly-ionised, 
metal-rich ejecta of the star, has a very high X-ray col 
umn density, typically about a few x 10 22 c m" 2 (Kovama et al. 



1994; Skinner et al. 2010] lZhekov et al.l 1201 U iGosset et al 



201 It lOskinova et al.ll2012l) . that is dust-free (or at least very 
dust-poor compared to the typical ISM). The typical sizes of 
WR nebulae (a few pc) are consistent with the maximum sizes 
estimated for the io nised GRB absorber (IWatson et al.ll2007t 
ISchadv et al.ll2.Ql lT though under a slightly different set of as- 
sumptions). The wind is already quite hot before the burst and 
so the bulk of the absorbing material is already significantly 
ionised and may be far enough away from the star that we 
would rarely if ever expect to see substantial transient ionisa- 
tion by the GRB itself. 

We need not confine our considerations to mass-loss sce- 
narios as observed for known, well-observed Galactic WR 
stars. It appears from recent work that many evolved massive 
stars suffer significant mass loss in outbursts, the archetype 
of which is ^Carinae, where the mass loss from the 1840s 
event has been es timated to have been as much as 15 M 
(Smith et al. 2003). These outbursts come, presumably, from 
stars transitioning via precisely this mass- loss mechanism 
to the WR phase dSmith & Owockil 12006). The nebulae 
around these high-mass stars are often dusty and very mas- 
sive, and appear to be more common than prev iously believed 
dWachter et alj|2010t iGvaramadze et al.ll2Q10b . Furthermore, 
such extreme mass loss episodes have been discovered in scat- 
tered light, dust reheating and absorption observed via su- 
pernovae, suggesting that eruptive mass loss is c ommon in 
the last millenium right up to explosion of the star dFox et all 
2010t iPastorello et al.ll2007t iSmith et al.ll20TTt lChandra"etaLl 



20121) . 

The progenitor ejecta scenario explains many of the avail- 
able observations, including the high level of the X-ray ab- 



sorbing column density compared to A^hi and its relative lack 
of evolution with redshift. 

However, in spite of the case made for progenitor ejecta, a 
major problem remains with this interpretation: for distances 
of the absorber from the progenitor larger than ~ 1 pc, the 
mass required becomes very large, however the luminosity of 
the prompt GRB and early afterglow strips light metals en- 
tirely of their electrons at anything except large distances. 

This minimum distance is so large that for a given column 
density, the total mass req uired is excessive for a progenitor 
wind. Detailed analyses (Perna & Lazzatil 12002) as well as 
our own simple photon number calculations, indicate that the 
GRB emission strips light metals to distances of several tens 
of pc in some cases and a few pc even for lower luminosity 
bursts. But at distances of 1 pc, as we note above, the mass 
required to produce the mean observed GRB X-ray absorp- 
tion is more than 1M Q in metals for spherically distributed 
ejecta. This value rises as the square of the distance, so that at 
a mean distance of even a few pc, we would typically require 
approximately ten solar masses of light metals. 

These large masses essentially mean r < 1 pc for progenitor 
ejecta, while in direct contradiction, the power of a typical 
burst completely ionises metals that lie closer than a few pc, 
thus exclud ing an effective absorb er at distance less than a 
few pc (e.g. lPerna & Lazz ati 2002). Highly asymmetric mass 
loss, with the mass ejected along the direction of the GRB, i.e. 
along the rotation axis of the star, could mitigate this problem, 
but not enough and we still need the mass to lie at distances 
typically smaller than ~ 1 pc. 

5. PROPERTIES OF THE ABSORBER 

We look at the properties listed above and examine specifi- 
cally the compatibility of these models with 1) the high level 
and distribution of Nh x , 2) the evolving nature of the dust-to- 
A'hx ratio, 3) the correlation between X-ray column density 
and the gas, 4) the large typical distances to neutral gas, 5) the 
absence of moderately-ionised species, and 6) any reports of 
features related to metals in the X-ray absorber. 

5.1. H II region or progenitor nebula ? 

Both the progenitor ejecta and natal H II region work well 
to reproduce the observed features of the X-ray absorption. In 
particular, the correlation between the X-ray absorbing col- 
umn density and the gas column density strongly suggests 
that the magnitude of the X-ray column is driven largely by 
confinement by the host ISM, i.e. for a given mass, the col- 
umn density through it decreases as the square of the ra- 
dius. As observed for WR nebulae and nearby H II regions, 
the radius of the region is related to galacto-centric distance 
dHunt & Hirashital 120091) . e.g. by far the la rgest WR ejecta 
nebula is above the Galactic plane (W71, Stock & Barlow 
2010). For a given ejected mass, a physically smaller bub- 
ble will result in a higher column density, giving rise to the 
relation we observe in Figs. [2] and [3] 

However both hypotheses suffer from difficulties. In the 
former scenario there is a 'mass-distance problem': a tension 
between the minimum distance the metals must lie at in order 
not to be stripped of all their electrons by the GRB and the 
maximum distance allowable to keep the mass required to re- 
produce the observed column densities at a reasonable level. 
We explore possible solutions to the mass-distance problem in 
the appendix, but find no solution consistent with the observed 
data and therefore reject it. In the latter scenario, the X-ray ab- 
sorbing column density should drop with redshift as the mean 



metallicity of star-forming regions drop. A He-dominated H II 
region absorption would resolve this issue. We explore this 
scenario below. 

5.2. Properties of a helium- dominated H II region absorber 

We examine here more quantitatively the prospect that the 
host H II region could have the properties required to repro- 
duce the X-ray absorption. 

The reason He comes to dominate the X-ray spectra of 
GRB afterglows is because the stripping radii related to a 
GRB are, somewhat counterintuitively, smallest for H, then 
He and then the metals considerably further out. This is be- 
cause the GRB emission has a far harder spectrum than even 
ho t stars; in photon fl ux terms, F(E) oc E~ r , where (V) ~ 
1 dKaneko et alj 120061) . The number density of the atoms 
thus dominates the ionisation distances, and so those ele- 
ments with higher number densities have the smaller stripping 
radii. Therefore we would expect H and He to feature promi- 
nently in the afterglow absorption, especi ally in low metallic- 
ity environments where G RBs are found (Modiaz et al.ll2008t 
iGraham & Fruc hter 2012). However, in H II regions, the hy- 
drogen is of course pre-ionised by the stars, while the 54 eV 
photons required to strip He are extremely rare, making He 
the dominant absorber. Indeed, we would expect a priori that 
He should absorb GRB afterglows very strongly and it's worth 
asking the question: if anything else were responsible for the 
excess X-ray absorption, then why don't we see He absorp- 
tion? 

We can calculate very roughly the effective stripping radius 
of various elements making somewhat simplistic assumptions 
equating the total flux of available photons and the total cross- 
sections of the hydrogen- and helium-like ions (i.e. with only 
one or two electrons remaining) at these radii. We have done 
this for a representative, w ell-studied GRB, GRB 050401 at 
z = 2.9 (Watson et al. 2006) and show these values in Fig. [4] 
It is clear that the space where the ionisation state of He re- 
mains largely unaffected is very significant compared to the 
stripping radius for O (the stripping radius for Fe is simi- 
lar to O, and is further out for C and other metals), leaving 
He as the only absorber in the region between about 5pc 
and ~ 30 pc. This marks the effective size of the H II re- 
gion surrounding GRB 050401. It must be larger than 5pc 
for there to be any He to absorb, and must be smaller than 
~ 30 pc, so that the metals and dust do not start to contribute 
substantially to the optical/UV spectrum. This limit on the 
size and our measurment of the column density also allows 
us to determine that the density of the H II region must be 
in the range 10 3 - 10 4 cm~ 3 . These numbers are consistent 
with the densities found for nearby H II regions, though we 
should note that GRB 050401 was a very luminous burst. In 
Fig. [4] we plot the effective X-ray column density expected 
from He II for the de nsities and sizes o f H II regions in the 
sample presented by iHunt & Hiras hita (2009). The mean 
columns for the radio-determined densities and sizes are sur- 
prisingly close the mean columns found for GRB afterglows 
(logA^x ~ 21.7). The column densities determined from the 
HST observations are consistently lower and may be affected 
by the res olution of the HST imag ing and ground-based spec- 
troscopy (IHunt & Hirashitall2009l . 

We have also plotted the therma l sublimation radius for dus t 
according to the prescription of Wax man & Drainel (2000). 
However, we note that in that paper they assumed only pho- 
tons in the range 1-7.5 eV would contribute to the dust heat- 
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Radius (pc) 

Figure 4. Equivalent hydrogen column density (No* ) as a function of 
radius for a sample o f nearby extragalactic star-forming regions from 
Hunt & Hirashita 1 2009) including many blue compact dwarf galaxies de- 
rived from radio observations (circles) and HST imaging with ground-based 
spectroscopy (squares). The distribution of the radio-derived column den- 
sities, while containing possibly severe selection biases, is compatible with 
the mean observed Nu x in GRBs I Watson & Jakobsson 2012; Campana et al. 
120121) , indicating that the bulk of the X-ray absorption in GRBs could 
arise from material in its surrounding H II region. The column density for 
GRB 050401 is indicated by a dashed line. The white region indicates the 
distance to which GRB 050401 has enough soft X-ray photons to fully strip 
helium and metals so that no X-ray absorption would be observed from the 
H II region with a radius smaller than this distance in an afterglow. Helium 
is stripped to smaller distances than oxygen (or other metals). Dust will be 
destroyed by thermal sublimation to large radii. The survival distances of 
He, O and dust from GRB 050401 are indicated by red, blue and shaded re- 
gions respectively. There is a large region where only He retains its electrons, 
making it the likely dominant X-ray absorber for most GRBs. 

ing since molecular and atomic hydrogen would absorb the 
UV photons above this energy. We use the parameters for the 
GRB 050401 burst from the Swift-BKT automated analysis^, 
as well as all the photons up to 40 eV, since there is no H2 or 
H I. This more than doubles the effective thermal dust subli- 
mation radius. The dust sublimation radius is rather close to 
the O and Fe stripping radii, and dust is therefore unlikely to 
play a role in the absorption by the H II region, again, leav- 
ing only He as the expected signature. The outcome of these 
considerations aligns comfortably with the fact that the op- 
tical/UV extinction observed in GRB afterglows is consistent 
with the column density of low-ionisation metals observed far 
from the burst (Zafar et al. in prep.) and with the fact that only 
very small column densities of moderat ely ionised metals ar e 
found in the spectra of GRB afterglows (Schady et al.ll20lTh . 

While H Lya is commonly observed in the optical spec- 
tra of z > 1.8 GRB afterglows, absorption lines associated 
with He are much more difficult to detect. The only lines 
we could reasonably expect to observe in absorption are far 
into the UV, at or above the He Lyman series. For singly- 
ionised He, this effectively means that the absorption is never 
observable, since for z < 3, the continuum photons in the He II 
Lya region of the spectrum are absorbed by neutral hydrogen 
in the Galaxy, and for z > 3, the continuum is absorbed by 
the He Lya forest. The proportion of He I to He II depends 
of course on the hardness of the radiation field ionising the 
region. However, for neutral He, a space-based, UV spectro- 

6 |http : / / gen . gsf c ■ nasa ■ gov/ notices_s/113120/BA/| 
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graph might find the He I resonance line related to the X-ray 
absorpton for GRBs at 1 < z < 2.5, dependent on the propor- 
tion of He I present. 

6. IMPLICATIONS OF HELIUM ABSORPTION 

We have concluded above that the data on X-ray absorption 
from GRBs allow us to exclude every scenario proposed to 
explain it to date with the exception of He in the natal H II 
region. Specifically, 

a) the absence of a strong decrease in the mean X-ray ab- 
sorption with redshift strongly implies that the X-ray 
absorber is not primarily due to metals in the ISM of 
the host galaxy; 

b) the correlation we find between Nn x and Ay implies 
an association between the X-ray absorber and the host 
galaxy; 

c) the evolution of the Nn x /Ay ratio coupled with the ap- 
parent lack of evolution of this ratio for the optical met- 
als to dust ratio (Zafar et al. in prep.) suggests that the 
X-ray absorber is not due to the metals in the dust that 
cause the observed Ay; 

d) that the evolution of the dust column is related to 
the metallicity evolution, while the X-ray absorber is 
largely independent of this; 

f) this in turn suggests that the X-ray absorber is due to 
He in the natal H II region or progenitor ejecta metals; 

e) the ionising power of the GRB precludes metals from 
the progenitor ejecta; 

g) the correlation between the X-ray absorbing column 
density and the dust extinction and the atomic hydro- 
gen column density indicates that the host galaxy ISM 
has an influence on the X-ray column density. 

We are thus left with a scenario in which the X-ray absorp- 
tion is caused by He absorption in the natal H II region that is 
undetected at UV/optical wavelengths and that is confined by 
the host galaxy ISM. Every other scenario proposed so far is 
excluded. 

This conclusion leads to some predictions about the observ- 
able properties of the absorber. First, as mentioned above, 
there might be a detectable signature of neutral He resonance 
absorption for GRBs in the redshift range 1 < z < 2.5, de- 
pending on the temperature of the H II region. Second, since 
the absorption is dominated by He, the absorption observed 
in X-rays should be smooth, i.e. absorption edges or lines due 
to metals should be weak. The claimed detection of S and Ne 
absorption edges in the afterglow of GRB 090 6 1 8 seems prob- 
lemati c in this respect (Camp ana et al.ll20TTl) . iCampana et al.l 
(1201 ll) found in their analysis that S and Ne were signif- 
icantly overabundant (though without Si or O being simi- 
larly overabundant, which might have been expected on nu- 
cleosynthetic grounds), but that the general metal abundance 
found was low. However, it may be possible that the Galactic 
foreground was underestimated: a n abundance low in metals 
(particularly in O) was assumed (Wilms et al. 2000) for the 
Galactic foreground absorp tion, substant ially lower than typ- 
ical for Galactic sightlines ( Watson [201 ll) . At the same time, 
the Ne I edge at z = 0.54, the redshift of the GRB, is very 
close in energy to the Galactic O I edge. Clearly, accounting 



for the Galactic O I edge is crucial to an accurate determina- 
tion of the z = 0.54 N e I edge. Furthermore, observationally, 
Camp ana et al.l (1201 ll) found the S and Ne abundances to be 
strongly correlated, such that if the Ne abundance was low, the 
S abundance would also likely be low. Therefore we believe 
that the detection of absorption edges in GRB 090618 may 
need to be revisited. No other non-transient absorption fea- 
ture has been reported in GRB afterglows, which is perhaps 
noteworthy in its own right, (though we are unaware of any 
study explicitly searching for these features in low-redshift 
GRB afterglows). Currently therefore, the data appear to be 
consistent with a smooth X-ray absorption. 

Fits performed previously to the X-ray spectra of GRB 
afterglows with solar metallicity material may therefore be 
somewhat unreliable not just in the absolute value of the 
derived equivalent hydrogen column density, as is generally 
acknowledged, but also in shape, especially at low redshift 
where the O absorption edge is especially prominent. The 
very low redshift GRB 060218, which is str ongly absorbed, 
has an evolving thermal comp onent reported (Campa na et al.l 
2006b; iGhisellini etail 120071) . The fits to this thermal com- 
ponent may be significantly affected by a He-dominated ab- 
sorber instead of a solar metallicity absorber. 

In addition to thermal components in low redshift bursts, 
the conclusions drawn about the metallicity of highly ab- 
sorbed bursts, especially at high redshift, based on the X- 
ray absorption assuming the absorption to be dominated 
by metals, should also be revisited in light o f our results 
( Cam pana et aDl2007t IWatson & Laursenll201 ll) . 

This first successful interpretation of the X-ray absorption 
allows it to be used for diagnostic purposes. Based on the lu- 
minosity of the burst and the continued existence of neutral or 
singly-ionised He, we can set lower bounds to the radii of the 
H II region in which the burst resides, as well as upper limits 
based on the fact that the metals must be stripped by the burst. 
With such limits, we can then convert the observed column 
density to volume densities for these regions. Ultimately, it 
should be possible to calculate approximate distributions for 
the densities and sizes of individual H II regions where mas- 
sive stars form across a vast range of redshifts, currently even 
as far as GRB 050904 at z = 6.3. 

7. CONCLUSIONS 

We have shown that the X-ray absorption in GRBs is cor- 
related with the dust extinction and with the neutral hydrogen 
column and that the ratio of Nn x /Ay changes with redshift 
in a way similar to the metallicity redshift evolution of GRB 
host galaxies. This suggests that the magnitude of the X-ray 
absorbing column density in GRBs is correlated with the gas 
density in the host galaxy ISM, explaining the relationship 
found between high Nh x bursts and so-called 'dark bursts'. 
Using these new findings, we excluded all models relating the 
X-ray absorber to anything outside the host galaxy including 
the warm-hot IGM. On the grounds that the Nn x /Av changes 
with redshift, we concluded that only helium in the host H II 
region or metals ejected by the progenitor star could be the 
primary X-ray absorber. However, the ionising power of the 
GRB sets a minimum distance for metals to retain any elec- 
trons and hence be effective X-ray absorbers; this places a 
minimum mass on the metals required which is too large for 
ejecta from the progenitor, excluding this hypothesis. We thus 
concluded that helium in the host H II region causes most of 
the X-ray absorption observed in GRB afterglows. This con- 
clusion allowed us to set limits on the sizes and densities of 
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the typical H II regions in which X-ray absorbed GRBs ex- 
plode, > a pc and < a few tens of pc, with densities about 
10 3 — 10 4 cm" 3 , consistent with observations of sizes and elec- 
tron densities of H II regions in the Milky Way and in nearby 
galaxies. 
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APPENDIX 

POSSIBLE SOLUTIONS TO THE MASS-DISTANCE PROBLEM FOR PROGENITOR EJECTA 

The mass-distance problem could possibly be resolved if the afterglow doesn't see the same absorption as the prompt emission, 
i.e. the material absorbing the afterglow is somehow not flash-stripped by the GRB. One apparently plausible scenario is that the 
afterglow has a substantially larger opening angle than the prompt emission. In this case, while the prompt emission may strip 
metals to very large radii, it does so along a very narrow beam, while the majority of the afterglow emission travels through H- 
and He-like-ionised material much closer to the jet head. Another possible solution might be to clump the ejecta. Finally, we 
could conceal the metals in the solid phase. In this latter case, optical absorption lines would not be observed. However, we 
would then need to resolve the lack of observed extinction and the dust destroying power of the GRB, both of which might be 
resolved by placing most of the metals in large grains. We explore these possibilities below. 

Implications of different prompt/afterglow opening angles 

A prediction of the different opening angles of the prompt and afterglow emission is that the prompt phase should have much 
lower column densities on average than the afterglow. Whether this applies to later X-ray flares is unclear, as they may have the 
same opening angle as the jet. Some fraction of cases where Swift-XRT follow-up of a very long burst, or a burst with a precursor 
is rapid enough, should show absorption clearly rising as the burst transitions from the pr ompt to afterglow phases. So far, 
howe ver, this has not been observed. Indeed, the opposite has been reported for GRB 050904 (Ge ndre et al.l2007l ICampana et alJ 
l2007h . and in our data we do not see clear evidence for a negative difference between prompt and afterglow absorption, in direct 
conflict with the requirements of this hypothesis. It therefore seems unlikely that the GRB prompt and afterglow emission see 
very different columns of gas. 

Other solutions worth considering for the mass-distance problem are that the X-ray absorber is clumped into extremely high 
column density knots, or extremely high volume density knots. In the former case, the knots have high enough column density 
that even at distances of 0. 1 pc, the burst is not powerful enough to strip all the electrons. The covering fraction must be extremely 
close to unity, however, and for the highest column density situations, the absorber must be close to the Compton optical depth. 
The extremely large covering fractions seem unlikely, making this scenario interesting but low probability. In the latter case, very 
high volume densities could potentially allow recombination on a short timescale (e.g. ~ 1 s). Such high volume densities, even 
at 0.1 pc, imply an incredible relative thickness of the shell (~ 10~ 10 ), indicating that this scenario is essentially impossible. 

Large grain hypothesis 

X-rays remove the inner electrons from metals. Once this happens, the electrons do not recombine in an observable time, 
even for naked nuclei. This is the problem we are trying to overcome. By placing the metals in the solid phase, the affected ion 
can recombine rapidly, and the energy in the liberated electron can be dissipated and radiated away at longer wavelengths. This 
process is more effective in larger grains. The question is whether any substantial fraction of the dust can survive the intensity 
of the GRB at distances of a fraction of a pc, sinc e it has been shown theoretically that the extreme radiation field of a GRB 
is sufficient to destroy dust to very large distances (Wax man & Drai ne 2000; [Fruchter et al.ll200"Tl iPerna & Lazzati 2002]). We 
hypothesise that the metals responsible for the X-ray absorption could be primarily in large (~ 1 /im) dust grains. First, such 
grains are more robust to UV sublimation. The UV heating is approximately proportional to the surface area of the grain. Larger 
grains have a higher volume to surface ratio and are thus more robust to sublimation. Se cond, such large g rains could solve the 
problem that we do not observe substantial reddening associated with the X-ray absorber (Pe rna et al.l l2003). However, on closer 
inspection, it appears they cannot. The absorption in X-rays has the characteristic spectral shape of neutral medium absorption. 
Specifically, the spectrum is not flat at low energies, indicating the covering fraction is close to unity. This means that if metals 
in large grains are responsible for the X-ray absorption, the grains must have a covering fraction close to unity. For most dust 
materials, the absorption efficiency is close to unity. Even for highly-transparent materials with relatively low absorption in the 
optical (e.g. diamond), and large grains, the scattering efficiency is close to unity, which would lead to very strong extinction 
of the optical and UV light. We thus conclude that we cannot find a way to reconcile the progenitor ejecta hypothesis with the 
observed properties of GRB afterglows. 



